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Abstract 

We study the hadronic D meson decays into a pseudoscalar meson P and an even-parity meson 
M, where M represents a scalar meson S, an axial-vector meson A, or a tensor meson T. These 
decays are first analyzed in the flavor-diagram approach. Fits to the SP modes with S being a 
non-strange scalar meson show that neither the simple qq picture nor the q 2 q 2 scheme is favored 
by data. Current measurements on the AP decays are insufficient for a meaningful analysis. 
Some TP data are inconsistent with the others. In certain cases, the H^-annihilation diagrams 
indicated by the data are unexpectedly large. As a comparison, we also compute their decay rates 
in the factorization approach using form factors extracted from the covariant light-front model. 
We find that factorization works well for Cabibbo-allowed D + — > SP, AP decays free of the weak 
annihilation contributions (lU-exchange or TU-annihilation) . For the other SP and AP modes, it is 
necessary to include weak annihilation contributions to account for the data. However, factoriztion 
fails for D — > TP decays for some unknown reason; the predicted rates are in general too small by 
at least two orders of magnitude compared to experiment. We also examine the finite width effects 
of resonances. Some decay modes which are kinematically forbidden become physically allowed due 
to the finite width of the resonance. We show that the branching fraction of D + — > an + extracted 
from three-body decays is enhanced by a factor of 2, whereas B(D° — > j 2(1270) K°) is reduced by 
a factor of 4 by finite width effects. 
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I. INTRODUCTION 



A plethora of interesting but puzzling phenomena regarding strong interactions have been re- 
vealed in the study of charmed meson decays. The magnitudes and phases of CKM factors associ- 
ated with the dominant decay amplitudes in such processes are well-determined. Moreover, these 
decays do not receive significant corrections from penguin-type loop diagrams. Such salient features 
enable one to readily extract the magnitudes and relative strong phases of various flavor diagrams. 
More specifically, there are four types of flavor digrams that dominate the charmed meson decays. 
They are the color-allowed tree amplitude, the color-suppressed tree amplitude, the T^-exchange 
diagram, and the VF-annihilation diagram. Previous studies [iHsl have shown the importance of the 
VF-exchange and ^-annihilation diagrams, presumably due to significant final-state rescattering 
effects, in the two-body decays of charmed mesons to two pseudoscalar (PP) or one vector and one 
pseudoscalar mesons (VP). Moreover, these amplitudes are seen to have nontrivial relative strong 
phases, also a result of sizable final-state interactions (FSIs). Among various sources of FSIs, the 
most important one is arguably the contribution from intermediate resonance states near the D 
meson masses. 

In this paper, we set to study the hardonic decays of charmed mesons into a pseudoscalar meson 
P and an even-parity meson M . Here M can be a scalar meson, denoted by S, an axial-vector 
meson, denoted by A or a tensor meson, denoted by T. The D — > SP decays have been studied 
previousl y in Refs. [9Ml5l]; the D — > AP decays in Refs. and the D — > TP decays in 

Refs. j25l427l| . In these decays, the flavor diagram of each topology has two possibilities: one with 



the spectator quark in the charmed meson going to the pseudoscalar meson in the final state; and 
the other with the spectator quark ending up in the even-parity meson M. We thus need two 
copies of each topological diagram to describe the decay processes. Many of these decays have 
been observed in recent years through dedicated experiments and powerful Dalitz plot analysis 
of multi-body decays. An extraction of the sizes and relative strong phases of these amplitudes 
therefore becomes possible. 

One purpose of studying these decays is to check our understanding in the structures and 
properties of light even-parity mesons. Another goal is to learn the FSI pattern in view of the rich 
resonance spectrum around the D meson mass range. 

This paper is organized as follows. In Section [XTJ, we review current experimental status of the 
measurements of multi-body charmed meson decays that are relevant to our analysis. We provide 
the information of flavor SU(3) classification, decay constants, and form factors for the light S, 
A, and T mesons in Section IIII1 Section IIVI presents the so-called flavor-diagram approach to the 
decays. Each decay mode is decomposed in terms of quark diagrams characterized by their flavor 
topologies. Current experimental data are used to infer the magnitude and strong phase associated 
with each of the amplitudes as best as we can. Under the factorization assumption, we compute 
the rate of each decay mode in Section [Vj We also examine the finite width effects for certain decay 
modes in Section IVIl Some discussions on the form factor model used in this work, theoretical 
uncertainties, comparison with other approaches are presented in Section [VIII A summary of our 
findings is given in Section [Villi 
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II. EXPERIMENTAL STATUS 



It is known that three- and four-body decays of heavy mesons provide a rich laboratory for 
studying the intermediate state resonances. The Dalitz plot analysis is a very useful technique 
for this purpose. We are interested in D — > MP decays extracted from the three-body decays of 
charmed mesons. Many results are available from ARGUS, Femilab, CLEO, FOCUS and BaBar. 
The results of various experiments are summarized in Tables Hl lIIII where the products of B{D — > 
MP) and B(M — > P1P2) are listed. To extract the branching fraction for D — > MP, we apply the 
narrow width approximation 

T(D -> MP PiP 2 P) = r(D -> MP)B(M PiP 2 ) . (1) 

(Finite width effects in certain decays will be discussed in Section EB) For the branching fractions 
of two-body decays of even-parity mesons, we shall use 

jB(/ (980) -> vr+vr") = 0.35 ± 0.08 , B(a+(980) -> r/vr) = 0.845 ± 0.017 , 

0(A^°(143O) -»■ A+vr") = ^(0.93 ± 0.10) , £(A^+(1430) -»■ A+vr°) = -(0.93 ± 0.10) , 

B(/ 2 (1270) -> vrvr) = (84.8±f;f)% , B(/ 2 (1270) -> M) = (4.6 ± 0.4)% , 

£(a 2 (1320) -»• M) = (4.9 ± 0.8)% , B(JfJ(1430) Avr) = (49.9 ± 1.2)% , (2) 

and B(/o(1500) -> vr+vr-) = §(34.9 ± 2.3)%, where we have applied the value of r(/ (980) -> 
7T7r)/[r(/ (980) -»• vrvr) + r(/ (980) -»• A A)] = 0.52 ±0.12 obtained by BaBar H and the Particle 



Data Group (PDG) average, r(a (980) -> AA)/r(a (980) -> vrr/) = 0.183 ±0.024 [29[]. To obtain 
the branching fraction of /o(980) —> vrvr, we have assumed that vrvr and A A are the dominant decay 
modes of f (980). 

Several remarks are in order: 

1. There are two measurements of — > /o(980)vr+ from —> vr+vr+vr~ by E687 |30J] and 



B(Dt -> /o(980)vr+)£(/ -> vr+vr") = \\ \ , ^ 2 (3) 



E791 [31| with the results 

' (6.3 ±0.8) x 10- 3 E791, 
(1.2 ±0.2) x 10" 2 E687. 

These two data are not used by PDG for the average. The most recent Dalitz plot analysis 
of D+ -»• vr+vr+vr- by BaBar yields B(D+ -> (vr+vr- ) s _ W ave7r+) = (0.92 ± 0.07) x 10~ 2 0. 
Since the 5-wave is the sum over /o(980), /o(1370) and /o(1500), it is clear that the E687 
result is too large. Hence, we will only quote the E791 result in TableUfor L>+ — > /o(980)vr + 



with /o(980) — > vr+vr . For an early theoretical study, see 33]. 

2. Many of the 3-body decays listed in Table U involve the decays /o(980) — > K + K~ and 
ao(980) — > K + K~ . Since the central values of the /o(980) and ao(980) masses are below the 
threshold for decay into a pair of charged kaons, the narrow width approximation Eq. (pTJ) is 
no longer applicable because /o or ao need to be off-shell. That is, the relation 

B(D° -> / (980)vr°; /o -> K+K~) = B(D° -> / (980)vr°)B(/ -> K+K~) , (4) 

does not hold. From Table HI we see that B(D° — > /ovr°;/o — > vr+vr- ) is smaller than 
B(D° — > /ovr°;/o — > K + K~) by one order of magnitude. Since B(fo — > K + K~) is smaller 
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than B(fo — > tt + tt~) due to phase space suppression 1 , the branching fraction B(/o(980) — > 
K + K~) obtained from B(D° — > /ovr°; /o — > K + K~) under the narrow width approximation 
will be too large by at least one order of magnitude. 

The decay D + — > K~ir + ir + is dominated by the S'-wave (K~ir + ) component which consists 
of -?Tq(800) (or k), ^(1430) and non-resonant contributions. The PDG value is B(D + — > 



(K ir + )s-w!ive'K + ) = (7.62 ±0.25)% 29]. However, PDG does not take the measurements of 
D + ^°(800)vr + and ^°(1430)vr + by E791, E691 and E687 for the average. If we take 
the E791 results alone [3lJ, we find 

B{D+ -> Ao*°(800)vr+) = (6.5 ± 1.9)% , B(D+ -> K*°(U30)ir + ) = (1.8 ± 0.3)% . (5) 

Hence, the decay D + — > ^q°(800)7t + has the largest branching fraction among the two-body 
D — > SP decays. 

Since the axial- vector mesons decay into three pseudoscalar mesons via strong interactions, 
their resonant substructures are studied in the Dalitz plot analysis of four-body decays. For 
example, information on the decay rates of D° — > i^ _ a^(1260) and -KT-f (1270)-7r + can be 
extracted from the study of D° — > K~ 7r + 7r~7r + . 



Because the a meson is very broad in its width, of 0(600 — 1000) MeV [29J, the use of the 
narrow width approximation is not justified and it becomes necessary to take into account 
the finite width effect of a. We will examine the finite width effect for the decay D — > air in 
Section EO 



III. PHYSICAL PROPERTIES OF EVEN-PARITY MESONS 
A. Scalar mesons 

It is known that the underl ying structure of scalar mesons is not well established theoretically 
(for a review, see e.g. Refs. 0, 38]). Many scalar mesons with masses lower than 2 GeV have been 
observed, and they can be classified into two nonets: one nonet with masses below or close to 1 
GeV, namely, the isoscalars /o(600) (or a), /o(980), the isodoublet ^(800) (or k) and the isovector 
ao(980); and the other nonet with masses above 1 GeV, namely, /o(1370), ao(1450), i^o(1430) and 
/o(1500)// (1710) 2 . If the scalar meson states below or near 1 GeV are identified as a conventional 
low-lying + qq nonet, then the nonet states above 1 GeV could be excited qq states. 

In the naive quark model, the flavor wave functions of the light scalars read 

a = —=(uu + dd) , /o = ss , 

V2 



1 The ratio B(f ->• K+K-)/B(f tt+tt") was measured to be 0.69 ± 0.32 by BaBar 28] and 0.25±g " 
by BES Q. 

2 Since not all three isosinglet scalars /o(1710), /o(1500), /o(1370) can be accommodated in the qq nonet 
picture, it is widely believed that one of them should be primarily a scalar glueball. 
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TABLE I: Experimental branching fractions of various 
venience, we have dropped the mass identification for 
taken from Ref. 29] unless specified otherwise. 



D — > SP decays. For simplicity and con- 
/ (980), a (980) and #5(1430). Data are 



B(D -> SP) x B{S -> P 1 P 2 ) 



B(D -> SP) 



B(D+ -s- f n+) = (4.5 ± 1.4) x 10- 4 

B(D+ -> / (1500)tt+) = (4.7 ±1.7) x 10" 4 

B(D+ -> f K+) = (1.6 ± 1.0) x 10- 4 
B(D+ -> <ttt+) = (2.1 ± 0.2) x 10~ 3 
B(D+ -> ft°JC+) = (1.0±g"|) x 10- 3 



B(D+ -> / 7r+)B(/o -> 7r+7r-) = (1.55 ± 0.33) x 10~ 4 
B(D+ -> / o (137O)7r+)0(/o(137O) -> tt+tt~) = (8 ±4) x 10" 5 
B(D+ -> /o(1500)7r+)B(/ (1500) -* Tr+Tr") = (1.1 ± 0.4) X 10~ 4 
B(D+ -* / (1710)tt+)B(/o(1710) -> tt+tt") < 5 X 10~ 5 
B(D+ -> f K+)B(f -> tt+tt ) = (5.6 ± 3.4) x lO" 5 
B(D+ -> o-tt+)B(o- -> TT+TT") = (1.37 ± 0.12) X 10" 3 
B(D+ -> k°K+)B(k° -> X-TT+) = (6.8t 3 ;^) x 10" 4 
B{D+ -> K* °K+;K* ° -> E-TT+) = (1.83 ± 0.35) x 10" 3 
B(D+ -> (X-TT+) s _ wavc TT+) = (7.62 ± 0.25)% 



W-) = (2.76l°;S)xlO- 3 



• /oK )B(/o 

• /o(1370)¥°)B(/ (1370) -> tt+tt") = (5.0 ± 1.2) x 10" 
► / 7T )B(/o -> TT+TT ) = (3.6 ± 0.8) x 10" 5 
' /ott°; fo -> X+if") = (3.5 ± 0.6) x 10~ 4 
. / (1370)tt )B(/o -> tt+tt ) = (5.3 ± 2.1) x 10" 
. /o(1500)tt°)B(/o 



B(.D 
B(D° 
Bp 
B(D° 
B(D° 
B(D° 
B(D° 
B(D° 
B(D" 
B(D° 
B(D° 
B(D° 
B(D° 
B(D° 
B(D° 
B(D° 
B(D° 
B(D° 
B(D° 
B(D° 



■ W+IT ) = (5.6 ± 1.6) X 10 

) = (4.5 ± 1.5) x 10" 



— > TT 7T" 



. /o(1710)tt°)B(/ 

► foK°; f -> K+K-) < 2.0 X 10~ 4 
. fo(1370) K°)B(f -> K+K~) = (3.6 ± 2.2) x 10" 



B(£>° -> f IC) = (8.0ll;|) x 10~ 3 
E(£>° -> /ott°) = (1.0 ± 0.3) x 10- 4 



B(£>° -> /o(1500)tt°) = (2.4 ± 0.7) x 10~ 4 



'' a 



• a+(1450)E" 
>%K+;a,o 



■ K+K°) = (1.24 ± 0.36) x 10 

. K+K°) = (1.5 ± 0.1) x 10 
. — n 



3 a 



K+K°) = (2.1 ± 0.1) x 10" 

1 



K-~K°) < 2.4 x 10 



agK°;og -> if+Jf - ) = (6.2 ± 0.8) x 10" 
a°K°;a° -> K+K~) = (5.3 ± 0.4) X 10" 

• agK K) )B(a° -> »;tt ) = (1.34 ± 0.42) X 10" 2 

• a-K°)B{a -» tt+tt") = (1.18 ± 0.21) X 10~ 4 

. k;-tt+)b(k;- 



3 

3 a 
12 



■K° 



— *o 



,+0.80) x 1Q - 



) = (4.90_ a64/ 
_ft:-7r ) = (4.6 ±2.1) x 10 
K* \ )B(Kq' -» i^-TT+) = (5.7±f-g) x IP' 3 



B(Z>° -> ag"^) = (1.6 ± 0.5)% 
B(D° -> <ttt°) = (1.8 ± 0.3) x 10" 4 
B(D° -> X*-tt+) = (7.91J-J) x 10~ 3 
B(£)° -> K*-n+) = (1.5 ± 0.7)% 
B(D° ^^"tt ) = (gjtl^j x IP' 3 



B(Ds 
B(Ds 
B(Dg 
B(Dg 
B(Dg 
B{Dt 
B{Dj 

B{Dt 



(TT+TT- ) S -wavc7T+) = (0.92 ± 0.07) X lO" 2 6 
/ TT+)B(/ -> TT+TT") = (6.3 ± 0.8) X 10" 3 c 

/ott+;/o -> Jf+Jf - ) = (6.0 ± 2.4) x 10" 3 
/ott+; /o -> Jf+Jf-) = (1.55 ± 0.13) x 10~ 2 d 

■ /o(1370)tt+)B(/ -> X+_ft:-) = (2.37 ± 0.35) x 10" 2 d 
/o(1710)tt+)B(/ -5- K+K~) = (1.87 ± 0.29) x 10" 2 d 
K* ° K+)B(K* ° -> _^-tt+) = (5.1 ± 2.5) x 10" 3 

■ ~K* ° K+)B(K* ° -> ^ _ tt+) = (2.15 ± 0.40) x 10" 3 d 
i^*°TT+)B(i<: * -5- K+-K-) = (5 ±4) X 10" 4 



B(D+ -> / 7T+) = (1.8 ± 0.5)% 



B(D+ -> K* °K+) = (8.2 ± 4.1) x 10" 3 
B(D+ -> ~K*oK+) = (3.5 ± 0.7) x 10" 3 
B(D+ -> Ar*°7T+) = (8.1 ± 6.5) x 10- 4 



"From 
b From 
c From 
d From 



BaBar 35 [. 
BaBas- 
E791 [an 
CLEO [3 



2]- 



Oq = -^(tm — dd) , a.Q = ud , a Q = du , (6) 

v2 

/t + = -us , k° = ds , R° = sd , kT = su , 

where the ideal mixing for fo and a is assumed as /o(980) is the heaviest one and a the lightest 
one in the light scalar nonet. However, this simple picture encounters several serious problems: 
(i) It is impossible to understand the mass degeneracy betweem /q(980) and ao(980). A related 
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TABLE II: Experimental branching fractions of D — > AP decays taken from Ref. 291 ] . 



Decay 




Experiment 


Decay 


Experiment 


D+ -> 


if?(1270)7r+ 


< 7 x 10~ 3 


D+ -> if°a^(1260) 


(7.0 ± 1.2)% 


D+ -> 


^(1400)tt+ 


(3.8 ± 1.3)% 


L>° -> if- a 7/(1260) 


(7.9 ±1.1)% 


£>° — > 


iff (1270)tt+ 


(1.14 ±0.32)% 


D° -»• if°a?(1260) 


< 1.9% 


L> ^ 


iff (1400)7r + 


< 1.2% 


L>° -»• 7r-<+(1260) 


(8.98 ±0.62) x 10~ 3 


L> ^ 


^f(1400)7T° 


< 3.7% 







B(D° -»• if+(1270)if T )£ (iff (1270) -»• if+vr+vr-) = (8.1 ± 1.8) x 10~ 4 
B(D° -»• iff ( 1400) if+ jiff (1400) -> if+vr+Tr") = (5.4 ± 1.2) x 10~ 4 



TABLE III: Experimental branching fractions of various D — > TP decays. For simplicity and 
convenience, we have dropped the mass identification for /2(1270), d2(1320) and if | (1430). Data 
are taken from Ref. {29] ] unless specified otherwise. 



B(D -» 


TP) x B(T -> PiP 2 ) 




B(D - 


+ TP) 




B(D+ - 


-)• / 2 VT + )e(/2 -)■ 7T+7r _ ) = 


(5.0 ±0.9) x 10" 4 


B{D+ 


-> /2VT+) 


= (8.8 ± 1.6) x 10" 4 


B(D + 


-> lfJV + )B(if? ) -> if-7r- 


f ) = (2.1 ±0.4) x 10~ 4 


B(D+ 




) = (6.3 ± 1.2) x 10~ 4 


B{D + 


-> if£°7r+)£(if|° -> if+TT 


~) = (5.0 ±3.4) x 10~ 5 


B{D+ 


-> ifl°7T 4 


•) = (1.5 ± 1.0) x 10~ 4 


B{D+ 




= (i.7±53) x 10- 4 




-> 4^°: 










B(D+ 


< 3.0 x 10- 3 



B{D° -> M°)£(/ 2 -> vr+vr-) = (1.91 ± 0.20) x 10" 4 £(D° -»• / 2 7r°) = (3.4 ± 0.4) x 10" 4 
B(D° -> / 2 if°)£(/ 2 -> vr+vr-) = (2.8±f;§) x 10" 4 £(D° -»• / 2 if°) = (5.0+|f) x 10" 4 

s(£>° -> k;-tt+)b(k;- -> ifV) = (7.ol 4 ;°) x 10- 4 b(d° -> if 2 *-^ + ) = (2.ii^) x 10- 3 



fi(L>° -> a+if-) < 2 x 10" 3 
B(D+ f 2 7r + )B(f 2 TT+7T-) = (1.1 ± 0.2) x 10~ 3 a £(£>+ ->■ / 2 vr+) = (1.9 ± 0.4) x lO" 3 
B{D+ -> if|°^+)S(if|° -> if+71--) = (5 ± 4) x 10~ 4 £(£>+ -»■ if|°7r+) = (1.5 ± 1.2) x 10~ 3 

a From BaBar [3J|. 



question is why ao is heavier than k if it does not contain a strange quark? This is the so-called 
"inverted spectrum problem." (ii) The P-wave 0+ meson has a unit of orbital angular momentum 
which costs energy around 500 MeV. Hence, it should have a higher mass above rather than below 
1 GeV. (iii) It is hard to explain why a and k are much broader than /o(980) and ao(980). (iv ) 
The 77 widths of oo(980) and /o(980) are much smaller than naively expected for a qq state 39]. 
(v) The radiative decay <p ao (980)7, which cannot proceed if ao(980) is a pure qq state, can be 
nicely described in the kaon loop mechanism [40j]. This suggests a considerable admixture of the 
KK component. 

It turns out that these difficulties can be readily resolved in the tetraquark scenario where the 
four-quark flavor wave functions of light scalar mesons are symbolically given by [4l[ 

a = uudd , /o = ss(uu + dd) / v2 , 
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ctg = —=(uu — dd)ss , ciq = udss , a = duss , 
V2 

k + = usdd , k° = dsuu , R° = sduu , kT = sudd . (7) 

The four quarks q 2 q 2 can form an 5-wave (not P-wave!) + meson without introducing a unit 
of orbital angular momentum. Moreover, color and spin dependent interactions favor a flavor 
nonet configuration with attraction between the qq and qq pairs. Therefore, the + q 2 q 2 nonet 
has a mass near or below 1 GeV. This four-quark description explains naturally the inverted mass 
spectrum of the light nonet, especially the mass degeneracy of /o(980) and ao(980). The fall-apart 
strong decays a — > tttt, k — > Kir and /o, «o — > KK are OZI super-allowed without the need of any 
gluon exchange. This explains the broad widths of a and k, while /o(980) and ao(980) are narrow 
because of the suppressed phase space for their decays to the kaon pairs. The decays of /o(980) and 
ao(980) are dominated by /o(980) — > tttt and ao(980) — > r/ir, respectively. Lattice calculations have 
confirmed that ao(1450) and Kq(1430) are qq mesons, and suggested that a and k are tetraquark 



mcsoma 



431 ] . Since exotic 4-quark states have not been seen experimentally, this may imply 
the structure of diquark-antidiquark bound states for the light scalar mesons (for a review, see 
Ref. |3). 

In the 2-quark picture with ideal mixing, /o(980) is purely an ss state. This is supported by the 
data of Df — > /ovr + and <p — > /07, implying the copious /o(980) production via its ss component. 
However, there also exists some experimental evidence indicating that /o(980) is not a pure ss state. 
First, the observation of T(J/ip — > /quj) ~ ^(J/ip — > /o<^>) 29] clearly shows the existence of the 



non-strange and strange quark contents in /o(980). Second, the facts that /o(980) and ao(980) have 
similar widths and that the fo width is dominated by tttt also suggest the composition of uu and 
c?J pairs in /o(980); that is, /o(980) — > tttt should not be OZI suppressed relative to ao(980) — > ttt]. 
Therefore, isoscalars a(600) and /o(980) should have a mixing 

|/o(980)) = |ss)cos6» + |nn)sin6» , |cr(600)) = -\ss) sm9 + \nn) cos 6* , (8) 

with nn = (uu + dd)/y/2. Experimental implications for the fo-a mixing angle have been discussed 
in detail in Ref. [13]: the mixing angle lies in the ranges of 25° < < 40° and 140° < 6 < 165° 3 . 

Likewise, in the four-quark scenario for light scalar mesons, one can also define a similar fo-a 
mixing angle 

|/o(980)) = \nnss) cos cp + | uudd) sin (p , |cr(600)) = — | nfiss) sin <\> + \uudd) cos <p . (9) 
It has been shown that (j) = 174.6° 47]. 



3 Recently CLEO has measured the semileptonic decay £)+ — > /o(980)e + ^ e with the result B(Df — > 
/o(980)e+^ e )S(/o -> 7T+7T-) = (0.20 ±0.03 ±0.01)% [45]. Using the value B(/ n+ir-) = (50^)% 
inferred from the BES measurement and the QCD sum rule prediction B(Df — > foe + is) = 
cos 2 6 x (0.41)% [46|, CLEO then extracted the mixing angle to be cos 2 9 = 0.98to;°i. However, this 
is subject to two major uncertainties. First, the branching fraction of /o(980) — > ir + ir~ has not been 
measured directly. For £>(/ (980) tt + it~) ps 0.35 as used in this work, B(Df / (980)e+^ e ) and cos 2 8 
would be enhanced by a factor of 1.4. Second, the theoretical prediction of this semileptonic decay is 
model dependent as it depends on the form factor of the Df — > /q transition with /q = ss. 
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In principle, the 2-quark and 4-quark descriptions of the light scalars can be discriminated in 
the semileptonic charm decays. For example, the ratio 

B(D+ -> f £ + v) + B(D + -> a£+v) 



R 



B{D- 



(10) 



is equal to 1 in the 2-quark scenario and 3 in the 4-quark model under the flavor SU(3) symmetry 
48J. In reality, the light scalar mesons may have both 2-quark and 4-quark components. Indeed, 
a real hadron in the QCD language should be described by a set of Fock states each of which has 
the same quantum number as the hadron. For example, 



a+(98o)> = rjud) + va>^> + CL><M + 



udg 1 



udss 



(11) 



In the tetraquark model, _ 3> ' l p a % while it is the other way around in the 2-quark model. 
The decay constant of the scalar meson is defined as 4 



(S(p) | <?27m<?1 1 0) = fsP, 



I-' ! 



(S|<? 2 <?i|0) = m s fs 



(12) 



The neutral scalar mesons a, /o and a[j cannot be produced via the vector current owing to charge 
conjugation invariance or conservation of vector current: 

U = ff = f a = ■ (13) 
Applying the equation of motion to Eq. (fl~2j) yields 

Vsfs = Is , with [i s 



m s 



m 2 (//) - mi{^) 

where mi and m\ are the running current quark masses. Therefore, the vector decay constant of the 
scalar meson f$ vanishes in the SU(3) or isospin limit. The vector decay constants of i^g(1430) and 
the charged ao(980) are non- vanishing, but they are suppressed due to the small mass difference 
between the constituent s and u quarks and between d and u quarks, respectively. The scalar 
decay constants fs have been computed in Ref. 49J] within the framework of QCD sum rules. 
From Eq. ^ we obtain / ao(980) ± = 1.0 MeV, / ao(145 o)± = 5.3 MeV, and /^ (143 o) = 35.9 MeV. In 
short, the vector decay constants of scalar mesons are either zero or small. 
Form factors for D — >• P,S transitions are defined by 5C ] 



(Pip'W^Dip)) 



P„ 



m 



D 



m. 



% ) pr^) + 



77? 



D 



7777 



(S(p')\A„\D(p)) 



where P.. 




77? 



D 



777t 



r 



9m 



pDS 



(Q 2 ) + 



in 



D 



(15) 



fj, — {p+p')fii Qfi = (p~p')tM- As shown in Ref. [51(, a factor of (— i) is needed in the D — ^ S 
transition in order for the D — >• S form factors to be positive. This can also be checked from heavy 
quark symmetry consideration 511 ] . 

Throughout this paper, we use the 3-parameter parametrization 

F(0) 



1 



a(q 2 /m 2 D ) + b(q 2 /m? D ) 2 



(16) 



for D — > M transitions. The parameters F(0), a and b for D — > S transitions calculated in the 
covariant light-front (CLF) quark model are exhibited in Table HVl 



For pseudoscalar mesons, the decay constant is defined as (P(p)|g27/j 75<7i |0) = —ifpp^- 



S 



TABLE IV: Parameters in the form factors of D,D S -> / (980), Kq( 1430) transitions in the 
parametrization of Eq. (116p . as obtained by fitting to the covariant light-front model 511. The 
numbers in parentheses are the form factors at q 2 = obtained using the ISGW2 model 521 ] . 



F F(0) a 6 


F F(0) a b 


F(f /09 0.49 (0.13) 0.07 -0.03 
F Dir*(i430) 48 ^ Q ^ _ Q n Q2 


F® sfoa 0.46 (0.23) -0.29 0.07 



B. Axial-vector mesons 

In the quark model, two nonets of J p = 1 + axial-vector mesons are expected as the orbital 
excitation of the qq system. In terms of the spectroscopic notation 2S+1 Lj, there are two types of 
P-wave axial- vector mesons, namely, 3 Pi and 1 P\. These two nonets have distinctive C quantum 
numbers for the corresponding neutral mesons, C = + and C = — , respectively. Experimentally, 
the J PC = 1++ nonet consists of ai(1260), /i(1285), /i(1420) and K XA , while the l+~ nonet 
contains &i(1235), /ii(1170), /ii(1380) and K\b- The physical mass eigenstates ifi(1270) and 
iTi(1400) are mixtures of the K\ A and K\b states (we follow PDG 29J to denote the 3 Pi and 1 P\ 
states of K\ by K\ A and K\b , respectively) , 

Ki(1270) = K lA sm9 Kl + K 1B cos8 Kl , 

#i(1400) = K 1A cos6 Kl - K 1B sm6 Kl . (17) 

Since these states are not charge conjugation eigenstates, consequently, mixing is not prohibited. 
Indeed, the mixing is governed by the mass difference between the strange and non-strange light 
quarks. There exist several estimations on the mixing angle 9k 1 in the literature. From the 
early experimental information on masses and the partial rates of iTi(1270) and -ftTi(1400), Suzuki 
found two possible solutions, each with a two-fold ambiguity, \9ki_ \ ~ 33° and 57° [53J]. A similar 
constraint 35° < ^ 55° was obtained in Ref. 5J] based solely on two parameters: the mass 
difference between the a\ and b\ mesons and the ratio of the constituent quark masses. An analysis 
of t — > K\(\270)v T and Xi(1400)^ r decays also yielded the mixing angle to be ~ 37° or 58° with 
a two-fold ambiguity [55]. Most of these estimations were obtained by assuming a vanishing fx 1B • 
With the help of analytical expressions of /k 1a is obtained in the CLF quark model [5l|, two 
solutions for the i^i(1270)-i^i (1400) mixing angle, 50.8° and —44.8°, have been found in Ref. 561 ] . 
However, the second solution has been ruled out by the measurements of B — > K\ (1270)7 an d 
B —7- K\ (1400)7 [56l |. Therefore, we shall use Qk\ = 50.8° in the ensuing discussions. 



For the decay constants and the form factors of the axial vector mesons, we shall follow Ref. 51] 
to define them as 5 

(A(p,e)\A„\0) = f A m A el , 



5 The relative signs of the decay constants, form factors and mixing angles of the axial-vector mesons 
were often very confusing in the literature. As stressed in Ref. [57j ]. the sign of the mixing angle 6k 1 



is intimately related to the relative sign of the K\a and K\b states. In the CLF quark model [51| and 
in pQCD [58|, the decay constants of K\a and K\b are of opposite signs, while the D(B) — > K\a and 
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TABLE V: Parameters in the form factors of D — )■ ai(1260), b\ (1235), Kia, K\b transitions in the 
parametrization of Eq. (116j) . as obtained by fitting to the covariant light-front model fo"l||. The 
numbers in parentheses are the form factors at q 2 = obtained using the ISGW2 model [52]. As 
noticed in the footnote of this subsection, the form factors for D — t- 3 P\ and D — > l P\ transitions 
are of the same (opposite) signs in the CLF (ISGW) model. 



F 


no) 




a 


b 


F 


F(0) 


a 


b 


1/ Dai 
V 

v 


0.31 (- 
0.34 (- 


-0.60) 
-0.37) 


0.85 
1.44 


0.49 
0.15 


V Dbl 


0.49 (0.64) 
0.44 (0.50) 


0.89 
0.80 


0.28 
0.27 



(A(p,e)\Ap\D(p D )) 
{A{p,e)\V^\D{p D )) 



m D - tua 



PdP 



P A DA {q 2 ) 



i{ (m D - m A )e* fl V 1 DA (q z ) - (e* ■ p D )(p D +p). 



yDA ( „2 



m D - tua 



-2mA 



e -Pd 



(18) 



Because of the charge conjunction invariance, the decay constant of the 1 P\ non-strange neutral 
meson such as 65(1235) must be zero. In the isospin limit, the decay constant of the charged 
b\ vanishes due to the fact that b\ has an even G-parity and that the relevant weak axial-vector 
current is odd under G transformation. As for the strange axial vector mesons, it is known that 
the decay constant of the 1 P\ meson vanishes in the SU(3) limit [3]. 

In the following, we shall take f ai = 238 ± 10 MeV obtained using the QCD sum rule method 
611 ] . similar to the p meson, f„ ~ 216 MeV. This means that the ai(1260) meson can be regarded 



as the scalar partner of the p meson, as it should be. In the CLF quark model 51[, if we increase 
the constituent d quark mass by an amount of 5 db 2 MeV relative to the u quark mass, we find 
/& 1 = 0.6 ± 0.2 MeV for the charged &i which is very small. Using the experimental results 
B(t -»■ liri(1270)i/ T ) = (4.7 ± 1.1) x IO- 3 and T(r -> #i(1270)i/ r )/[r(r -»■ #i(1270)z^ T ) + T(r -»■ 



#l(1400)r/ T )] 



0.69 ± 0.15 

l/ifi(1270)l : 



291 ]. we obtain 



169.5±|f:| MeV 



1/^(1400)1 — 

In the CLF quark model the signs of the decay constants Jk 1a 
MeV and f KlB 



139.2l^J MeV 



12 MeV 



(19) 

and fx 1B are fixed: /k 1a = —212 
561 ]. This together with the mixing angle 6ki = 50.8° also fixes the 



signs of fx x to be 

M(1270) = "170 MeV , / Kl( i400) = "139 MeV , (20) 

where we just consider the central values. 

Finally, the D — > A form factor parameters in the CLF quark model are given in Table |V] 



D(B) — > K\b form factors are of the same sign. The mixing angle Bk x is positive. It is the other way 
around in the approaches of QCD sum rules [59( and the ISGW model [52], |60(: the decay constants of 
K\a and K\b have the same sign, while the D(B) — > Kia and D(B) — > K\b form factors are opposite in 
sign. These two conventions are related via a redefinition of the K\a or K\b state, i.e., K\a ~K\a or 
K 1B -> -K 1B . 
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C. Tensor mesons 



The observed J p = 2+ tensor mesons / 2 (1270), $(1525), a 2 (1320) and #2(1430) form an SU(3) 
1 3 P 2 nonet. The qq content for isodoublet and isovector tensor resonances are obvious. Just as the 
rj-rf mixing in the pseudoscalar case, the isoscalar tensor states f 2(1270) and /2(1525) also have a 
mixing, and their wave functions are defined by 

/ 2 (1270) = -L(/« + / 2 d )co S %+/|sin0 /2 , 

/a(1525) = -±=(f2 + fi) S me h -flco S e h , (21) 

with /| = qq. Since tttt is the dominant decay mode of ,f 2 (1270) whereas /2(1525) decays predomi- 
nantly into KK (see Ref. 29 H, it is obvious that this mixing angle should be small. More precisely, 



it is found that 6 h = 7.8° [62] and (9±1)° [29(. Therefore, / 2 (1270) is primarily an (uu + dd)/y/2 
state, while /2(1525) is dominantly ss. 

The polarization tensor e^ v of a 3 P 2 tensor meson with J PC = 2 ++ satisfies the relations 

e„ u = e vll , = , p^ v = Vv eT = , (22) 

where is the momentum of the tensor meson. Therefore, 

(0|(F - A) ll \T(e,p)} = ae^tf + be? vPlt = , (23) 

and hence the decay constant of the tensor meson vanishes identically; that is, the tensor meson 
cannot be produced from the V — A current. 

The general expression for the D — > T transition has the form 60] 

(T(e,p T )\(V - A) fl \D(p D )) = ih(q 2 )e^ pa E*» a p Da {p D +p T y(p D -p T y + % 2 )e>£ 

+ H(q 2 )^PDPD(PD+PT)» + b-(q 2 )elpp a D pP D (p D -p T )v. (24) 

The form factors h, k, b + and 6_ have been calculated in the ISGW quark model [60j and its 
improved version, the ISGW2 model 521 ] . They are also computed in the CLF quark model 51] 
and listed in Table IVII 

The decay amplitude of D — > TP always has the generic expression 

A(D TP) = E%py D M(D -> TP) . (25) 

The decay rate is given by 

T(D -> TP) = ( \M(D -»■ TP)\ 2 , (26) 

where p c is the magnitude of the 3-momentum of either final-state meson in the rest frame of the 
charmed meson. 



IV. DIAGRAMMATIC APPROACH 

It has been established sometime ago that a least model-dependent analysis of heavy meson 
decays can be carried out in the so-called topological diagram approach. In this diagrammatic 
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TABLE VI: Parameters in the form factors of 
parametrization of Eq. (116j) . as obtained by fittin. 
form factor k is dimensionless, while k, b + and 6_ 
theses are the form factors at q 2 = obtained usin 



D -> a 2 (132), #1(1430) transitions in the 
g to the covariant light-front model [5j|. The 
are in units of GeV -2 . The numbers in paren- 
g the ISGW2 model [52]. 



F 


F(0 




a 


b 


F 


F(0) 


a 


b 


h Da 2 


0.188 


(0.203) 


1.21 


1.09 


k Da 2 


0.340 (0.613) 


-0.07 


0.12 


b + 


n HQ A 

— U.Uo4 


v — U.UOZ J 


u.y i 


U.Oo 


L.Da 2 


n i on { a naA \ 
U.lzU (^11.1)04] 


1.10 


U.DD 


h Dh q 


0.17 


(0.20) 


1.28 


0.90 


k Df 2q 


0.27 (0.61) 


-0.21 


0.12 


b D hq 


-0.08 


(-0.05) 


1.02 


0.51 


b D hq 


0.10 (0.06) 


1.14 


0.49 


h DK 2 


0.192 


(0.14) 


1.17 


0.99 


k DK* 


0.368 (0.71) 


-0.04 


0.11 


K 


-0.096 


(-0.060) 


1.05 


0.58 


b DK5 


0.137 (0.069) 


1.17 


0.69 


h Dshs 


0.15 


(0.21) 


1.04 


0.79 


k D s f 2s 


0.59 (1.15) 


0.22 


0.09 


b+ shs 


-0.09 


(-0.09) 


0.95 


0.54 


b D s f 2s 


0.13 (0.12) 


1.05 


0.60 


h D s K* 2 


0.15 


(0.35) 


1.11 


0.99 


k D s K* 


0.42 (1.16) 


0.08 


0.11 


b + 


-0.07 


(-0.08) 


0.96 


0.60 


h D s K* 


0.11 (0.13) 


1.10 


0.63 



scenario, all two-body nonleptonic weak decays of heavy mesons can be expressed in terms of six 
distinct quark diagrams 63H65I]: T, the color-allowed external VF-emission tree diagram; C, the 
color-suppressed internal VF-emission diagram; E, the ^-exchange diagram; A, the VK-annihilation 
diagram; P, the horizontal VF-loop diagram; and V, the vertical I^-loop diagram. (The one-gluon 
exchange approximation of the P graph is the so-called "penguin diagram".) It should be stressed 
that these diagrams are classified according to the topologies of weak interactions with all strong 
interaction effects encoded, and hence they are not Feynman graphs. All quark graphs used in this 
approach are topological and meant to have all the strong interactions included, i.e., gluon lines 
are included implicitly in all possible ways. Therefore, analyses of topological graphs can provide 
information on FSIs. Various topological amplitudes in two-body hadronic D decays have been 
extracted from the data in 0, [s, 66-72] after making some reasonable approximations, e.g., flavor 
SU(3) symmetry. 



The topological amplitudes for D — > SP, AP,TP decays have been discussed in [13|, l22j, |27| ]. 
There are several new features. First, one generally has two sets of distinct external IF-emission and 
internal VF-emission diagrams, depending on whether the emitted particle is an even-party meson 
or an odd-parity one. Let us denote the primed amplitudes T' and C for the case when the emitted 
meson is an even-parity one. Second, because of the smallness of the decay constants of even-parity 
mesons except for the 3 Pi axial-vector state, it is expected that |T'| <C |T| and |C"| <C \C\. This 
feature can be tested experimentally. Third, since Kq and the light scalars a, k, /o, ao fall into 
two different SU(3) flavor nonets, in principle one cannot apply SU(3) symmetry to relate the 
topological amplitudes in D + — > /ovr + to, for example, those in D + — > Kq°tt + . 
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TABLE VII: Topological amplitudes and branching fractions for various D — > SP decays. In 
Scheme I, light scalar mesons a, k, ao(980) and /o(980) are described by the qq states, while Kq as 
excited qq states. In Scheme II, light scalars are tetraquark states, while Kq are ground-state qq. 



The /o — o mixing angle 9 in the 2-quark model is defined in Eq. ([8]). The experimental branching 



fractions for D° 


-> #o~7r+ and ^ ^o°^ + 


are 


taken from Table U after averaj 


r P For siinnliritv 


we do not consider the /o — <J mixing in the 4-quark model. 




Decay 


Amplitude (I) 




Amplitude (II) 


*3expt 


D° -> f n° 


^V c * d V ud (-C + C'-E-E') sine 




±V* d V ud (-C + C -E-E') 


(1.0 ±0.3) x 10" 4 


— > f K 


V 2 




+v* a v ua c 




V* 'K w f^(C + £)sin<9 + E' cos01 




A=V*V„ r ?(C + 2E' + E) 


(8.0±%l) x 10" 3 


— r CLq A 






V C s V ud\l- +£') 










V* 3 V ud (C - E)/V2 


(1.6 ± 0.5)% 


-> a~K+ 






V* d V ua (T + E') 




% _ 
— > aj7r 


V^ d (T' + i?) 




v; d v ud {T' + E) 




— > 0,q1T + 






V*V ud (T + E') 




— > air 


|V;* d Vnd(-C + C — E — E') cos 




±V; d V ud (-C + C>-E-E>) 


(1.8 ±0.3) x 10" 4 










D+ -+ / 7T+ 


^V c * d K d (T + C' + A + A') sine 

+v* a v us c' cose 




■^V c * d V ud (T + C + A + A') 

+V2V C * v us c 


(4.5 ± 1.4) x 10" 4 




V* d V ua (T + A') sin 6 + A cos 6] 




^V;* d V us (T + 2A + A') 


(1.6 ± 1.0) x 10- 4 


->-a°7r+ 


V? s Vud(T' + (7) 

V c * d V„ d (-T -C-A + A')/V2 




V*V ud {T' + C) 

V c * d V ud (-T -C'-A + A')/V2 






-^V c * d V ud (T + C + A + A') cos 6 
-V* S V US C sine 




v; d v ud (T + c + A + A') 


(2.1 ± 0.2) x 10" 3 




V c * s V U3 T + V* d V ud A 




V£V u3 T + V c * d V ud A 


(i.oi°:«) x io- 3 




V* s V ud {T cos e + (A + A') sin 

V c * V us [(T + C" + A) cos + -±=A' sin 




VZV ud (2T + A + A')/V2 


(1.8 ±0.5)% 


-> /oi^+ 




■j=V* a V us (2T + 2C + 2A + A') 

+^v; d v ud c 






V* s V ud (T + E') 




V*V ud (T + E') 


(8.2 ± 1.4) x 10" 3 




V c * s V ud (C - E')/V2 




V* 3 V ud (C'-E>)/V2 


(9.2^-1) x IO" 3 




V c * s V ud (T + C) 




V* s V ud (T + C) 


(1.8 ±0.3)%° 


D+ -> K* °K+ 


V c * s V ud (C + A) 




V c * s V ud (C'+A) 


(3.6 ±0.7) x IO" 3 




v; d v ud T+v ca v* s A 




v; d v ud T + v cs v* a A 


(8.1 ± 6.5) x 10" 4 



a Data from E791 [3l|; see also Eq. ©. 



A. D^SP 

The topological amplitudes for D — > SP decays are listed in Table IVTTl for two different schemes. 
In Scheme I, light scalar mesons a, k, ao(980) and /o(980) are described by the ground-state qq 
states, while Kq as excited qq states. In Scheme II, light scalars are tetraquark states, while Kq 
are ground-state qq. The expressions of topological amplitudes are the same in both Schemes I 
and II except for the channels involving /o or a. 

Since the decay constant of fo and a vanishes, one can set T' = C = 0. From Table IVIII we 
have 

B(D + -» /Q7T+) = B(D° -> /qtt ) _ f tan 2 9 2-quark, 

B(D+^a7r+) B(D°^air°) \ I 4-quark. 1 ' 

It appears that the data of D + — > /o7r + ,cj7r + favor the 2-quark picture of the light scalars, while 
the measurements of D° — > /o7r°,o"7r° prefer the 4-quark scenario. Moreover, the D° — > clqK~ 
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TABLE VIII: Extracted flavor amplitude parameters from fits to the D — > SP decays, where S 
only refers to lighter scalar mesons here. The amplitude magnitudes are in units of 10 -6 GeV. In 
these fits, we set C = 0, E' = E, and A' = —A, with reasons explained in the text. The strong 
phases 5e,a are associated with the E and A amplitudes, respectively. For Fits (A) and (B), we 



take the mixing angle 


9 = 25°. 










Scheme 
Fit 


(A) 


I 


(B) 


II 

(A) 


(B) 


\T\ 
\A\ 
Sa 


2-i4±8;lI 

o i«+o.io 
(31 ± 2)° 






1.55 ±0.07 

o i c;+0.30 
z - io -0.45 

(351?!)° 




\C\ 
\E\ 
5e 






i 2n +0 - 05 

± - zu -0.06 

(168^)° 




1 qn+0-36 
i.yu_ 22 

i 1O+0.05 
i - i °-0.06 

(152jJ) 


Xmin / d - oi 


2.61/2 




5.07/1 


8.29/2 


0.74/1 



and ag 7r modes will be dominated by the IV-exchange diagram, E. The D+ a^K mode 
is dominated by the C amplitude. The D° -> K* °tt°, D + K^7r+, and D+ K$°K + are 
dominated by the E' , T, and A amplitudes, respectively. 

Table IVIII is divided into two parts separated by double lines. The upper part involves only 
light scalar mesons (/o, ao, <r, and k), whereas the lower part involves the Kq mesons in the heavier 
nonet representation. This division is made because the amplitudes of the same topology in these 
two groups have no a priori relations. We first note that none of the currently measured modes 
involve the T" amplitude. Secondly, one can simplify the upper part of the table by setting C = 0, 
for the decay constants of scalar mesons are expected to be either identically zero or relatively 
small. Moreover, the modes in the lower part of Table WTH have the same amplitude decomposition 
in the two schemes and involve the T, C, E' , and A amplitudes. One cannot set C = here 
because the decay constant of Kq is non-negligible, as commented after Eq. (|14|) . In this case, 
there are more theory parameters than observables, barring a fit. 

In the following, we will perform two sets of fit [(A) and (B)] to the modes involving only the 
lighter scalar mesons in the flavor diagram formalism. Fits (A) includes the five measured D + and 
Df decays. As noted in Eq. (|27p . the /ovr + and <T7r + modes are related because of the same flavor 
amplitude combination. Here we have to assume a relation between A and A 1 in order to reduce 
the number of parameters. Without further theoretical guidance, we have tried the cases A' = A 
and A' = —A for simplicity and found that the latter renders an equally good or better fit than the 
former. The strong phase 5a in Table IVIIII is measured with respect to T, which is assumed real. 
We also note here that all the strong phases given in the table are subject to a two- fold ambiguity 
(5 — > —5). From the Xmin va hies, one sees that Scheme I fits better than Scheme II in these modes. 
This is understandable because the 2-quark picture explains better the observed rates of the /ott + 
and an + decays. In either scheme, |^4| is about 1.5 times larger than \T\, showing the importance 
of the TV-annihilation contribution. Besides, the extracted relative strong phase is robust. 
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TABLE IX: Topological amplitudes and branching fractions of D — > AP decays. The notation is 
explained in the main text. Theory predictions are made within the factorization approach, with 
the mixing angle 9k 1 = 50.8°. 



Decay 








Amplitude 








Theory 


Experiment 


D+ -> 


if?(1270)7r+ 




1(Ta 


+ C' A )smO Kl + {T B 


+ c B ) 


cos 




4.7 x 10~ 3 


< 7 x 10~ 3 


D+ -> 


^?(1400)tt+ 






+ C' A )cos9 Kl -{T B 


+ C B ) 


sin 


o Kl ] 


2.2% 


(3.8 ± 1.3)% 


D+ -> 


K°af (1260) 






V*V ud (T' A + C A ) 








8.2% 


(7.0 ± 1.2)% 


D+ -> 


i?°6^(1235) 






V£V ud (T' B + C B ) 








2.2 x 1(T 3 




D° — >■ 


iff (1270)vr+ 


VcsVud 


[(Ta 


+ E' A )sm9 Kl + (T B 


+ E' B ) 


cos 




5.2 x 10~ 3 


(1.14 ±0.32)% 


D° — >■ 


iff (1400)vr+ 


v c * s v ud 


[{Ta 


+ E A )cos6 Kl - (T B 


+ E' B ) 


sin 




1.4 x 1(T 4 


< 1.2% 



- - 1 \ OS t * tJ> L \ -** 1 — J\J ^ 1 

D° -> X?(1270)tt° ^Kd^K^-^sin^ „ £ , 

-C -> X?(14fl0)7T° VSKrf^ [(C A - cos 0^ - (C' B - E' B ) sin 9 Kl ] 3.2 x 10" 3 < 3.7% 

D° if+(1270)if- [(T^ + £ A ) sin Kl + (3* + E B ) cos 9 Kl ] 4.6 x 10~ 4 

L>° iff (1270)if + F MS [(T A + E^) sin Xl + (T B + £ B ) cos 9 Kl ] 8.2 x 10~ 5 

£>° if 1 ± (1270)if =F 5.4 x 10~ 4 (8.1 ± 1.8) x 10" 4 
£>o if- a +(i260) y c * Kd(r A + E A ) 2.7% (7.9 ±1.1)^ 

J D°^/f°a?( 



if°a?(1260) V c * s V ud ^(C A - E A ) 1.2 x 10~ 4 < 1.9% 

n-af (1260) V;* d Kd(^ + E A ) 5.1 x 10~ 3 (8.98 ± 0.62) x 10~ 3 



£>o ^ ^-6+(1235) V*V ud (T B + E B ) 1.7 x 10" 

£>o ^ ifo 6 o (1235) O^^Cb - E B ) 3.0 x lO" 4 



Fit (B) includes the four measured D° decays. Here it does not matter what relation we assume 
between E and E', as far as the Xmm va l ue is concerned. The only effect is on the size and phase of 
the E amplitude. This is because the fair and air modes are related by the same flavor amplitude 
combination, as also noted in Eq. (|27p . We assume E' = E for an explicit fit. In Table IVITlj the 
strong phase 5e is measured with respect to C, which is assumed real. The Xmin values show 
that Scheme II explains this set of data better, as noted below Eq. ()27p . It is worth noting that 
the magnitudes and relative phase extracted in either scheme are roughly the same. Also, the 
amplitudes C and E are almost opposite in phase, as required primarily by the D° — > a^K decay. 
Finally, based on the Xmin vames °f these fits, the current data still cannot differentiate the two 
schemes yet. 



B. D -> AP 

The topological amplitudes for D — > AP decays are given in Table IIXI Instead of using sub- 
scripts of A and P to complicate the notation, we use the primed (unprimed) amplitudes to indicate 
that the spectator quark in the D meson ends up in the pseudoscalar (axial-vector) meson in the 
final state. The subscripts A and B refer to the amplitudes associated with the 3 Pi and x Pi 
axial- vector mesons, respectively. However, as we will see later in Section IV Bl the factorization 
approach predicts that such a distinction is only necessary for the T amplitudes. The assump- 
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TABLE X: Topological amplitudes and branching fractions of D — > TP decays. The notation is 
explained in the main text. Theory predictions are made within the factorization approach, with 
the mixing angle 9f 2 = 7.8°. 



Decay 






Amplitude 


Theory 


Experiment 


D + -> 


r -4- 


4^ V Md cos % (T + C" + A + A') 


0.9 x 10~ 6 


(8.8 ± 1.6) x 10 4 




4--rr0 

4k 




+V* s V us sme f2 C 






D+ -> 




V* s V ud (T' + C) 


4.5 x 10~ 7 


< 3.0 x 10~ 3 


D + —?■ 


K 2 7T + 




V* s V ud {T + C) 


1.9 x 10~ 5 


(6.3 ± 1.2) x 10 


D+ -> 






v; d v us (C + A) 





(1.5 ± 1.0) x 10~ 4 


D° — » 


/27T 


\V^V ud cos 9 h {C'-C-E'- E) 


5.1 x 10~ 8 


(3.4 ±0.4) x 10" 4 








+-^v*v usS me h c 






D ^ 




V*s V ud 


A= cos 9 h {C + E) + sin 9 h E' 


1.5 x 10~ 7 


(5.0±|;f) x 10- 4 


D° —?■ 






V* s V ud {T' + E) 





< 2 x 10~ 3 


D ^ 






V* s V ud {T + E') 


7.5 x 10" 6 


(2.1±J;?) x 10" 3 


Df^ 




v* s v ud 


A- cos 6 h (A + A) + sin 9 h T 


7.0 x 10~ 6 


(1.9 ±0.4) x 10~ 3 








V* cd V ud T + V* s V us A 


2.4 x 10~ 6 


(1.5 ± 1.2) x 10~ 3 



tion of Ca = Cb and Ea = Eb can be checked by comparing the rates of D° — > K°a5(1260) and 
D° — > K°bi(1235), which are seen to be roughly the same up to a tiny phase space difference. With 
the flavor symmetry assumption, the magnitudes of the invariant amplitudes of D° —> 7r~a^(1260) 
and D° — > K^af (1260) should differ by a factor of A ~ 0.2253, which is to be compared with 
0.171 ± 0.013 given by the current data. A distinctive feature between the Cabibbo-allowed D° 
and D + decays is that the H^-exchange diagrams (E) only involve in the former. 

Current data for D — > AP decays (only six branching fractions) are still insufficient for a 
sensible fit. The theory predictions in Table ILXl are based on the factorization calculations given 
in Section IV Bl 

C. D^TP 

The topological amplitudes for D — > TP decays are given in Table I2D Here we also use the 
unprimed (primed) symbols to indicate that the spectator quark of the D meson ends up in the 
even-parity (tensor) and odd-parity (pseudoscalar) mesons in the final state, respectively. There 
should be no confusion even though the amplitude symbols used here are identical to those in the 
D —> SP case. 

As described before, the decay constants of the tensor mesons vanish identically. Therefore, one 
can set T' = C = in Table [Xl The decay D + — > ' K + is kinematically forbidden as the K 2 
mass is above the kinematic threshold, though it is physically allowed through the width of K 2 . 
We therefore will not include it in our fit. Note that the measured D + and Df decays only involve 
the T, A, and A' amplitudes. The results of fits to these decay modes are given in Table IXll Since 
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TABLE XI: Extracted flavor amplitude parameters from fits to the D — > TP decays. The amplitude 
magnitudes are in units of 10 _6 GeV _1 , the same as M in Eq. (|26p . In these fits, we set C = 
and 8f 2 = 7.8°. The strong phases 5a,A' are associated with the A and A' amplitudes, respectively, 
relative to T. The contents of different fits are described in the main text. 



Parameter 


\T\ 


\A\ 5 A 


\A'\ 


8a> 


Xmin / d -°- f - 


Fit (A) 
Fit (B) 


i oe+0.14 
l-8t)_0.16 

O 1 O+0.39 

°- io -0.42 


12.63^;^ (0±3)° 
1.59 ±0.08 (173l?|)° 


11 24 +0 ' 16 


(180 ± 3)° 


11.4/0 
0/0 



there are only three measured modes in the D° decays whereas at least four additional parameters 
have to be introduced even if we set E' = E, it is impossible to determine the magnitudes and 
strong phases of C and E by considering a global fit to the TP decays. For that, a determination 
of B(D° — > a^K~) is crucial. 

Fit (A) in Table IXj includes all the available D + and Df decay modes except for D + — > K2*°K + ; 
there are thus 5 observables for 5 parameters. According to the results of Fit (A), the current 
data favor relatively large VK-annihilation diagrams. This is because \A\ and \A'\ are constrained 
respectively by the doubly Cabibbo-suppressed D + — > K^ir^ mode and the singly Cabibbo- 
suppressed Df — > mode to be large. Moreover, A and A' are about the same size but 

opposite in phase so that the branching fraction of Df —> /2vr + falls in the ball park. On the 
other hand, \T\ is largely constrained by the Cabibbo- favored D + — > K2*°tt + mode to be small. 
The largest contribution in the Xmin va hie comes from the D + — > /2vr + mode (~ 10.1). This is 
a manifestation of the disparity between the Cabibbo-favored Df — > /2vr + decay and the singly 
Cabibbo-suppressed D + — > f2^ + decay that are seen to have similar branching ratios, if \T\ is 
constrained not to play a role here. 

In view of the possibly problematic D + — > i^|°7r + and K 2 *°7r + modes, we exclude them in Fit 
(B) and set A' = A for simplicity; there are then 3 observables for 3 parameters. We note that it 
is not illuminating to consider A' = —A here because two of these modes involve the combination 
A ± A' as the major contribution. In Fit (B), \T\ becomes larger and A much suppressed. Also, 
the relative strong phase between A and T is almost opposite to that in Fit (A). 

The third column in Table |X] lists theory predictions based on factorization assumption to be 
discussed in Section IV CI A comparison between the predictions and the measured values shows 
an apparent deficit in theory account of the decay amplitudes. First, the magnitude of the tree 
contribution in the factorizationa approach is even smaller than the value of |T| in Table IXII 
Secondly, it is necessary to invoke the annihilation type of amplitudes to explain the observed 
data. 



V. FACTORIZATION APPROACH 



The diagrammatic approach has been applied quite successfully to hadronic decays of charmed 
mesons into PP and VP final states j^jls, 67-72]- When generalized to the decay modes involving 
an even-parity light meson in the final state, it appears that the current data are still insufficient 
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for us to fully extract the information of all amplitudes. Moreover, as shown in Tables IVIIII and 
IXIl the extracted parameters do not present a coherent picture yet. Therefore, we take the naive 
factorization formalism as a complementary approach to estimate the rates of these decay modes. 
In this framework, the IF-exchange and -annihilation type of cotributions will be neglected. We 
discuss the three categories of decays in the following subsections separately. 



A. D ^ SP 

The factorizable amplitudes for the D —> SP decays involve 

x (ds,p) = { p iq) \ {V _ A) ^ 0){S (p)\(V - Ar\D( PD )), 

X ( DP ^ = (S(q)\(V - AU0)(P(p)\(V - Ay\D(p D )), (28) 

with the expressions 

= _ fp{m 2 D _ m |)i^V) , * (DP ' S) = fs(ml - m%)F 1 DP (q 2 ) , (29) 
where use of Eqs. f)12|) and f)15j) has been made. The decay amplitudes of D — > KqP thus read 



A(D+ <V+) = ^ V ; s V ud [- ai U(ml-ml*)Fr o (<) 



+a 2 f K *{m' D -ml)F^{m 2 K , 
A(D°^K*-n+) = -^v; s V udai U(ml - m%*)F DK Z (ml) 



' 



A(D° ^ K* \°) = ^V? s V ud a 2 f K *(ml-ml)F D "(m 2 K *) , 
A(Df -+ K*°7r+) = ai V: d V ud U(m 2 Da - rn 2 K *)F°° K Hml) , (30) 

and likewise for the other D —> SP decays. 

Using the decay constants and form factors given in Section IIII Al the predicted rates of 
(D,D S ) — > (fo,KQ)P are computed and listed in Table IXIIl where we have used a\ = 1.22, 
0,2 = —0.66 and taken the form factors for D to ir and K transitions from the recent CLEO-c 



measurements of semileptonic D meson decays to tt and K mesons 73J. In order to test the 



factorization approach, we should focus on the modes in which weak annihilations (lU-exchange 
or VK-annihilation) are absent or suppressed. The Cabibbo-allowed decays D + — > Kq°it + and 
Df —7- /o 7r+ satisfy this criterion: the weak annihilation amplitude is absent in the former and 
suppressed by the fo — c mixing in the latter. We see from Table IXIIl that factorization works 
well for these two modes. For Cabibbo-allowed D + — > PP or VP decays, it is known that the 
color-allowed T and color-suppressed C amplitudes interfere destructively due to the opposite sign 
of the parameters a\ and a 2 - However, it is the other way around for Cabibbo-allowed D + — > SP 



and AP decays. From Eq. (|30p . it is obvious that the a± and a 2 terms in the decay ampli- 
tude of D + — > Kq°7t + interfere constructively. If they interfered destructively, one would have 
B(D + —J- K*qTt + ) = 5.4 x 10 -5 which is too small compared to experiment. Numerically, we obtain 
\T\ = 9.2 x 10~ 7 GeV and \C'\ = 1.0 x 10~ 6 GeV for this mode. Therefore, even though C is 
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TABLE XII: The predicted branching fractions for various D — > SP decays with the scalar mesons 
treated as qq ground states. For simplicity, we have dropped the mass identification for /o(980) 



and ifg(1430). The fo — cr mixing angle 6 is taken to be 25°. Theory predictions are made within 
the factorization approach in which the weak annihilation topologies (E and A) are neglected. 


Decay 




Amplitude 


^theory 


^expt 




-> 
D + -» 

Df -> 


/ovr + 


|V^K d (-C + C-E-E') sin 

KsKdl^ (C + sin + £' cos 0] 
■^V* d V ud (T + C + A + A') sm 9 

+V c * s V us C'cos6 
V *d V us[^(T + A') sin 6 + A cos 0] 
V* s V ud (T cos + (A + A') sin 0/>/2) 


7.8 x 1(T 6 

3.5 x 1(T 4 
1.4 x 1CT 4 

1.1 x 1(T 5 
1.3% 


(1.0 ±0.3) x 10" 
(8.0l||) x lO- 3 

/At" 1 ~l A \ t /~\ 

(4.5 ± 1.4) x 10 

(1.6 ± 1.0) x 10" 
(1.8 ±0.5)% 


-4 

_4 
-1 


D° — >■ 
-)■ 
D+ -> 


if 7T+ 

K *°vr+ 


V* s V ud (T + E>) 
V* s V ud (C - E')/V2 
V*V ud (T + C>) 
V* s V ud {C' + A) 
V: d V ud T + V cs V* s A 


2.0 x 1(T 3 
1.2 x 10 -3 
2.3% 

5.5 x 1(T 4 

2.6 x 1(T 4 


(8.2 ± 1.4) x 10" 
(9.2tH) x 10- 3 
(1.8 ±0.3)% a 
(3.6 ±0.7) x 10" 
(8.1 ±6.5) x 10" 


-3 

-3 
-4 



°Data from E791 [31|; see also Eq. ©. 



suppressed by the smallness of fx* (1430) an d a 2, it is enhanced sizably by the mass squared term 
(m 2 D — m 2 ) and the form factor at q 2 = m 2 K ». 

From Table IXLTl we see that the predicted rates for the other D — > KqP (P = ir, K) decays are 
smaller than experiments by a factor of 2 ~ 8. Note that they always receive weak annihilation 
contributions (E or A). Under the factorization hypothesis, the factorizable VF-exchange and 
^-annihilation amplitudes are suppressed due to the smallness of the form factor at large q 2 = 
m 2 D . This corresponds to the so-called helicity suppression. However, sizable long-distance weak 
annihilation can be induced via FSIs. For charm decays, it is expected that the long-distance weak 
annihilation is dominated by resonant FSIs. That is, the FSI via qq resonances is usually the most 
important one due to the fact that an abundant spectrum of resonances is known to exist at energies 
close to the masses of the charmed mesons. The diagrammatic-approach analysis in the last section 
suggests that weak annihilation diagrams are comparable to or even larger than the color-allowed 
tree amplitude T. Therefore, it is conceivable that the inclusion of E' and A amplitudes can 
account for B(D — > KqP). For D(Df) — > fair and JqK decays, the calculated branching fractions 
are typically too small by about one order of magnitude. To enhance the rates in this case, the 
weak annihilation contributions have to be larger than the color-allowed tree amplitude, as shown 
in Table IVIIIl Such an amplitude hierarchy poises a difficulty in theoretical understanding. While 
the non-strange content of /o(980) is small in the two-quark model for light scalars, it is not so in 
the tetraquark picture. This suggests that one should treat the light scalar mesons as bound states 
of qqqq. Unfortunately, the naive quark model is not applicable to evaluating the form factors for 
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the transition of D to a 4-quark state. 



B. D -> AP 



The factorizable amplitudes for the D — > AP decays involve 

x (da,p) = {PmV - A)tim A(p)\(V-Ar\D(p D )}, 

= (A( q )\(V-A),\0)(P(p)\(V-Ar\D( PD )). (31) 

with the expressions 

X( DA >V = 2f P m A V DA (q 2 )(e* ■ p D ) , = -2f A m A Ff p (q")(e* • p D ) . (32) 

It is then straightforward to write down the factorizable amplitudes of D — > i^i(1270)7r and 
D —7- i^i(1400)7r decays (dropping the overall e* ■ pz> terms for simplicity): 



A(D + -> ^(1270)vr+) = -^V c * s V ud 



A(D^ 



K°(1400)vr + ) 



2a2mK 1 (1270)fK 1 (1270)FF V ( m K 1 

^v c * s v ud 



2aim Kl ( 127 o)fn{sin6 Kl V 1A (ml) + cosO Kl V 

(1270) j 



'(ml)) 



2aim Kl(1Am) f n (cos6 Kl VQ KlA {ml) - sm6 Kl V DKlB (ml)) 



2a 2 m Kl ( 1400 ) f Kl ( 1400 ) ( m Ki 



(1400) i 



A(D° -> iff (1270)tt + ) 
yl(D° -> iff (1400)vr+) 
10 ' K°(1270)vr ) 



Gf t/* t/ 

^7f 



2aim Kl ( 1 270)/ 7 r(sin6'x 1 V ' DKlA (m^) + cos fl^ V^" 13 (ml)) 
2a l m Kl{ i A Q )f v (cosB Kl VQ KlA (ml) - sm6 Kl V DKlB (ml)) 



rDKi 



A(D 
A(D° 



Kl( 1400)tt° 



G F 



V* s Vud 2a 2 m Kl ( 1270 ) /jf 1 (1270) ( m /^ (1270) - 



(33) 



A(L>° -»• if+(1270)if~) 
A(L>° iff (1270)if + ) 



2 

~V2 



v*v 

v cs y us 



V2 cs us 



2a 2 m Kl ( 1400 ) f Kl ( 1400) F j 071 (mf^ (1400) ] 
2aim Kl ^270) /jd (1270) ^ ( m lr 1 (1270) 
2aim_ ftri ( 127 o)/^(sin6' i <r 1 F D 1A (m 2 K ) + cos 0^1^ 1B (m|-)) 



where we have taken into account the K\ A — if is mixing given by Eq. (|17p and neglected the 
short-distance factorizable VK-exchange contributions. Likewise, the D — >• If ai (1260) and D — > 
if b\ (1235) decay amplitudes read 



A(D + -+ IT af (1260)) 
A(D° K~at (1260)) 
A(D° -^if°a?(1260)) 
^(D° -> ir~ at (1260)) 
A(D + ^K°bt (1235)) 



2a 1 f ai m ai F 1 DK (m ai ) - 2a 2 f K m ai V uai (m' K 



Dai 1 



-^0^2ai/ 01 m 01 if *(f<) , 
-^V^V ud 2 ai f ai m ai F^(mlt) , 

-^v: s Vud r - - -- ^ 



2a 1 f bl m bl Fi J *(ml 1 ) - 2a2^m 6l Vo D6l (m^) 
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A ( D o ^ K-6+(i235)) = -^KXd 2a 1 f bl m bl F 1 DK (m 2 bl ) , 

A(D° ^K°6?(1235)) = ^^KA/™^™ 1 ^) • (34) 

Using the decay constants and form factors presented in Section IIIIBt the predicted rates 
of D — > Kiir, Ka±, Kbi decays are listed in Table HXl To test the validity of the factorization 
hypothesis, we focus on the D + decay to K\ir + and X°ai(1260) which are free of contamination 
from weak annihilations. We see that the predictions are in agreement with experiment for these 
D + decays and hence factorization works for D —> AP, just as the case of D — > SP. The 
predicted rates for D° decays are slightly smaller, which implies the importance of the TU-exchange 
contribution to D° decay modes. The theoretical calculations presented in Table IIXI are for the 
Kx(1270)-Ki (1400) mixing angle 6 Kl = 50.8°. When the other solution 6 Kl = -44.8° is used, 
we find the predictions B(D+ -> ^(1270)tt+) = 2.9%, B{D+ -> 1^(1400)^+) = 2.0 x 10~ 3 and 
B(D° — > K± (1270)ir + ) = 9.1 x 10 -5 , all in sharp disagreement with the data. Historically, it 
was first pointed out in Ref. pa!\ that a negative mixing angle Bk\ is ruled out by the data of 
D+ -> K?(1270)vr + and D° K{ (1270)tt + . 6 It was realized later that the negative Kl solution 
is also ruled out by the experimental measurements of B — )■ K\ (1270)7 an d B — > K\ (1400)7 UM- 

It is of interest to notice that the D° — > Kf : (14O0)K Zfl decay is not kinematically allowed, 
yet a branching fraction comparable to B(D° — > K^(1270)K T — > K T K ± ir~ir + ) has been ob- 
served. Since the width of Ki(1400) is 174 ± 13 MeV, the D° -> ^(1400)^ decay followed by 
1^(1400) -)• K ± tt + tt^ is certainly allowed. 



C. D^TP 



Since the decay constant of tensor meson vanishes, the factorizable amplitude of D 
always involves the expression 

X^ p ) = (P( q )\(V-A),\0)(T(p)\(V-Ar\D( PD )) 

= ifp^vPoPD k(mp) + b + (m 2 p){m 2 D - m T ) + b-(mp)mp 



TP 



(35) 



where use has been made of Eq. (|24p . The decay rate is given by Eq. (|26p . In general, TP final 
states are suppressed relative to PP states due to less available phase space. More precisely, 



T(D TP) 



2 p| / m D Y 

3 p p \TTIT J 



M(D TP) 



(36) 



r(D^PiP 2 ) 3p P \m T J M(D^PiP 2 ) 

where we have changed the notation slightly so that px denotes the c. m. momentum of the tensor 
meson and p p is the c. m. momentum of the pseudoscalar meson Pi or P 2 in the charmed meson 
rest frame. The kinematic factor \{pT/pp){ m D /^-t) 4 is typically of order (1 — 4) x 10~ 2 GeV~ 4 . 
An inspection of Table |X] indicates that, in the absence of weak annihilation contributions, the 
Cabibbo-allowed decays D + — > K* 2 G ^ + and D° — > should have the largest decay rates as 



6 The calculation of [55j was performed in the ISGW2 model [52] which has the opposite sign convention 
to the CLF model. 
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they proceed through the color-allowed tree diagram T. It is easily seen that all other TV-emission 
amplitudes in D — > C12K, D — > f 2 ir and D — > f 2 K are suppressed for various reasons. For example, 
it is suppressed by the vanishing decay constant of the tensor meson, by the small ji — f 2 mixing 
angle, by the parameter a 2 , or by the Cabibbo angle. 

From Table 1X1 we see that the predicted branching ratio of D + — > K* 2 ir + is of order 10 -5 , which 
is about two orders of magnitude smaller than experiment. Indeed, the theoretical calculation 
gives |T| = 3 x 1CP 7 in the unit of GeV -1 , which is much smaller than the value of \T\ listed 
in Table IXIl As for the decay D° — > K 2 ~ir + , its rate is similar to that of D + — > ^* 7r + but 
receives an additional W-exchange contribution. A fit of this mode to experiment will require 
\E\ > \T\, namely, VF-exchange dominates over the external IV-emission. The current measurement 
of B(D + — > K*2 ] K + ) is problematic as it is Cabibbo-suppressed and yet the measured rate is larger 
thanB(L>+ ^K* 2 \+). 

All the predictions shown in Table |X] are too small by at least two orders of magnitude, as 
originally noticed in Ref. [23]. In order to resolve the enormous discrepancy between theory and 
experiment for TP modes, one may consider possible form factor enhancement and finite width 
effects. One may compare the D — >• T transition form factors calculated in different models: the 
relativistic light-front quark model (see Table |V|) and the ISGW and ISGW2 quark models (see 
Table II of Ref. 23]). Since the form factors obtained from different models are of the same order, it 



is very unlikely that they can be enhanced by one order of magnitude to ameliorate the discrepancy. 
The finite width effect of the tensor resonances will be discussed in the next section. 



VI. FINITE WIDTH EFFECTS 

Normally we apply the narrow width approximation to extract the two-body branching fraction 
B(D — > MP) from 3-body decay data with M standing for an even-parity meson. There are three 
cases where the narrow width approximation is not valid or justified and the finite width of the 
resonance has to be taken into account: (i) The decay D — > MP is not kinematically allowed. For 
example, D + K* ° K+ ,K* 2 ° K+ and D° ( 1450) K~ , K^(im)K~ are forbidden if the scalar 
resonances are very narrow and on their mass shells, (ii) The resonance width is not negligible. For 
example, the widths of a and n are very broad, of order 600 — 1000 and 550 ± 34 MeV, respectively 



291 ] . (iii) The strong decay of resonance is marginally allowed or even forbidden kinematically. For 
instance, the central values of the /o(980) and oo(980) masses are below the threshold for decaying 
into a charged kaon pair. 



In general, the rate of the three-body decay D — > PiP 2 P is given by Ref. 75] 7 

2m D J( mi+m2 )2 27r 87rm| 



D 



1 2 AVV,m?X) 

g SP 1 P 2 v^'J 



(q 2 -m 2 s ) 2 + (T l2 {q 2 )m s ) 2 ^ FlF2 %W 



7 The case for D — > AP — > VP\P — > P1P2P3P is more complicated and has been discussed in Ref. 22[. 
The formula for T(D —> TP — s- P1P2P) given in Ref. [23] was erroneous and it is corrected here. 
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via a scalar resonance, and 



r(D _ TP _ w) _ i / ,m — >* £ | M(D _ rP) |. 4 



VT 7(mi+m 2 ) 2 27T 127T 

via a tensor resonance, where A is the usual triangular function A(a, 6, c) = a 2 + b 2 + c 2 — 
lab — 2bc — 2ca, mi (ro 2 ) is the mass of Pi (P 2 ), p(<? 2 ) = X 1 ^ 2 (m 2 J ,q 2 ,mp)/(2rnD), p'(q 2 ) = 
,mi,m2)/(2y/q 2 ~), and gup x P 2 is the strong coupling to be defined below. The "running" 
or "comoving" width Ty^a 2 ) is a function of the invariant mass mi2 = \/? °f the P1P2 system 
and has the expression 76] 

T 12 (g 2 ) = < mi2 V^"^ / "WSWW' ' (39) 

|r 5 ^4^T for M = S. 

^ ° m l2 p'(m 2 ) 

Note that the propagator of the resonance has been assumed to be of the Breit-Wigner form. From 
the measured widths of the resonances, one can determine their strong couplings 

r^p^l™^, r^ W = *ag»(£f. (4o> 

When the resonance width is narrow, the expression of the resonant decay rate can be simplified 
by applying the so-called narrow width approximation 

1 7T 

- m 2 M ) 2 + m 2 M T 2 M (q 2 ) " ^ " ' ^ 

It is easily seen that this leads to the factorization relation Eq. ([1]) for the resonant three-body 
decay. 

In the following, we illustrate the finite width effects with a few examples. 

D+ -» KpV+ -> JC + i^-7r + 

With a width of 270 ± 80 MeV for #5(1430), the decay D + followed by 

K + i^ _ 7r + is now physically allowed. In this case one should evaluate the two-step process T(D + — > 
K*qK + — > K + K~n + ) and compare the resonant three-body rate with experiment. Using Eq. (|37p 
and assuming that the coupling gsp 1 p 2 is insensitive to the variation in q 2 when the resonance is 
off its mass shell, we obtain 

B{D + -> K* °K+ -> K + K-ir + ) = (1.3±°;i) x 10~ 4 . (42) 

This is one order of magnitude smaller than the experimental value, (1.83 ± 0.35) x 10~ 3 (see 
Table U). Since 

A{D+ -> K* Q °K+) = V C * S V US T + V c * d V ud A , (43) 
the inclusion of VK-annihilation A will improve the discrepancy between theory and experiment. 

D + -> 
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Since the width of the a resonance T a = 600 — 1000 MeV is of the same order of magnitude as 
its mass, it is important to see its effect on the extraction of the branching fraction B{D + — > <77r + ). 
To see this, we shall first define a quantity 

_ T(D ^ MP ^ P X P 2 P) 
11 ~ r(D -> MP)B{M -> P 1 P 2 ) ' [ ' 

As r] goes to 1 in the narrow width approximation, the deviation of r] from unity gives a measure 
of violation in the factorization relation (TjQ). We first compute r\ theoretically; that is, both T(D — > 
MP — > P1P2P) and T(M — > P1P2) are computed in a model. The ratio 77 is independent of the 
form factor for D — > M transition. The factorization relation Eq. (TjQ) is then replaced by 

T(D -»■ MP -> PiP 2 P) = r? T(D -»• MP)B(M -> PiP 2 ) • (45) 

From the experimental input of /3(£> — > MP — >• P1P2P), we can then determine — > MP). 

For P + —7- cr7r + — >• 7r + 7r + 7r~ decays, we find 77 = 0.55 for r CT = 600 MeV and rj = 0.41 for 
T a = 1000 MeV, where we have taken B(a — > it + it~) = |. This means that the branching 
fraction B(D + —> air + ) = (2.1 ± 0.2) x 10 -3 listed in Table U should be enhanced by a factor of 
1/7/ = 1.8 ~ 2.4, depending on the a width. 

D z± TP 

We have examined the finite-width effect on the D — > TP channels listed in Table [Xj The 
measured decay widths of various tensor mesons are of order 100 MeV 
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The singly Cabibbo-suppressed decay D + — > K* 2 G K + — > K + K ir + is physically allowed due to 
the width T K ,o = 109 ± 5 MeV. From Eq. jHSJ) we obtain 

B(D + -> K* 2 °K + -> K + K~-k + ) = 4.3 x 10" 8 . (46) 



This is about 3 to 4 orders of magnitude below the experimental result (l.Tjlo'i) x 10~ 4 29]. 
However, this measurement seems to be problematic as the branching fraction of the Cabibbo- 
favored mode D + — > i^2°7r + — > K^ir + ir + is of the same order of magnitude, (2.1 ± 0.4) x 10~ 4 . 

For the other D — > TP decays, we compute the ratio 77 defined in Eq. (|44p and find that 
rj ~ 1.0 — 1.2 for most cases except for D° — > f^K^ where rj = 4.0. This means that the branching 
fraction of D° — > f2~K° extracted in Table ITTT1 should be reduced by a factor of 4 when the effect of 
finite width is taken into account 8 . 



VII. DISCUSSIONS 

A. The covariant light-front model 

We have relied heavily on the CLF model to obtain the form factors needed in this work. This 
relativistic quark model preserves the Lorentz covariance in the light-front framework and has 



been applied successfully to describe various properties of pseudoscalar and vector mesons [771 ] . 



The analysis of the CLF model has been generalized to even-parity, P-wave mesons in [5J|. Since 



8 Our conclusion for finite-width effects on D — >• TP differs from that in Ref. 



27|. 
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relativistic effects can manifest in heavy-to-light transitions at large to maximum recoil where 
the final-state meson becomes highly relativistic, the use of the non-relativistic quark model is 
probably not suitable here. Therefore, we believe that the CLF approach can provide more accurate 
behaviors of B — > M transitions at large recoil. For example, the tensor form factors for B — > K\ 
transitions at q 2 = derived in the CLF model lead to a prediction for B{B — > K\ (1270)7) m 
much better agreement with experiment than any other models (see Table IV of jo"^). 



B. Theoretical uncertainties 



Thus far we have not given the error bars to the theoretical results of decay rates. Here we 
discuss the possible sources of uncertainties in this study. In the CLF model we use the decay 
constants together with the given constituent quark masses to determine a fundamental parameter 
(3 in the model for describing the meson wave functions 5l|. For the decay constants of scalar 



mesons we use those obtained in 49J where QCD sum rules are employed and errors are given 
explicitly. For axial- vector mesons, the sum rule approach gives f ai = 238 ± 10 MeV 6l( while the 



experimental data of r — > K\v yield fK 1 (mo) = — (170±20) MeV and /^moo) = — (139±43) MeV 
[cf. Eqs. (fT9|) and ([20]) ]. The error bars will be propagated from decay constants to the parameter 
f3, and then finally to the form factors of interest. Therefore, the uncertainty analysis in the CLF 
model is quite involved and highly nontrivial. We will leave this task to a future publication. 



C. Comparison with other works 

In this work we have performed the study of the nonleptonic decays of charmed mesons in two 
approaches: flavor-diagram analysis and naive factorization. In the latter approach, form factors 
are obtained from the CLF quark model and the weak annihilation diagrams are neglected at the 
outset. In the literature most of the relevant studies are also based on the factorization approach, 
differing mainly in the values of the form factors to be used and the treatment of weak annihilation, 
n Tables IIVIIVII we have displayed in parentheses the form factors evaluated in the ISGW2 model 



521 ]. an improved version of the non-relativitsic quark model by Isgur, Scora, Grinstein, and Wise 
(ISGW) [60]. We see that the ISGW2 model predicts much smaller D — >• S form factors than the 
CLF model and other models not listed in Table ITVl In contrast, the form factors V® 0-1 , V® bl , 
k Da2 , k D ^, k Ds ^ and k D * K Z calculated by the ISGW2 model are much larger than the CLF 
results. 

As stressed before, the factorization hypothesis is best tested in the decays in which weak 
annihilation contributions are absent or suppressed. For the SP modes, the Cabibbo-allowed 
decays D + — > Kq°tt + and Df — > fo^ + belong to this category. It turns out that the former mode 
is ideal for testing different form factor models. The contribution to D + — > K 7r + from the color- 



suppressed tree amplitude C was not considered in [lOj, [ll|] , presumably due to the smallness of 
the Kq decay constant. Both contributions of color- allowed and color-suppressed amplitudes were 
taken into account in [9]. However, owing to the destructive interference in D + — > Kq°tt + and the 
smallness of 02, the predictions of T(D + — > K V) < T(D° K* 7T+) and the large suppression 
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of D° — > K*°ir° relative to D° — > K* ir + given in [9!] are not borne out by experiment. In the 
CLF model, the relative sign of the factorizable amplitudes X^ DS ' P ^ and X^ DP ' S ^ are fixed to be 
negative [see Eq. ([29]) ]. As a consequence, we conclude that, based on the CLF quark model, the 
interference in the Cabibbo-allowed D + — > SP decays must be constructive, contrary to the case 
of D -> PP, VP. 



Recently, the decays D —> K\n have been considered in 23|, [24| with D — > K\ from factors 



evaluated in the framework of the ISGW model and QCD sum rules, resepctively. While the 
predicted rates are similar to ours in most cases, the branching fraction of D + —> K®(l270)ir + 



was found to be (5.85 ± 0.37)% for 6 Kl = -37° and (3.18 ± 0.25)% for Kl = -58° in [24J. All of 
them are too large compared to the experimental limit 7 x 10~ 3 [29]. In our study we did include 



the color-suppressed tree amplitude [see Eq. ([33]) ]. which was considered in [23|] but neglected 



m 



24j . Since the decay constant of i^i(1270) is negative [Eq. (j20|) ]. it is clear from Eq. ([33 



that the ai term contributes destructively to D + — > i , T 1 (1270)'7r + . Consequently, our prediction 
B(D+ -)• K ( j ) (1270)7r + ) = 4.7 x 10~ 3 (Table EJ) is consistent with experiment. 9 

The above few examples indicate that the CLF model takes care of the relative signs of decay 
constants, magnitudes of form factors, and various hadronic matrix elements correctly. 



D. Comparison with B decays 

Charmful decays B — > D**P with D** = Dq , D\ , D[ , D\ denoting even-parity charmed mesons 
and charmless decays B — > MP (M being li ght even-parity mesons) have been extensively studied 



both experimentally [29] and theoretically 49j, l78l480l] . It is instructive to compare the present 
work with the B decays. 

While factorization works well for Cabibbo-allowed D + — > SP, AP decays, predictions are typi- 
cally about one order of magnitude smaller than experiment for the other decay modes, conceivably 
due to the negligence of weak annihilation contributions arising from final-state interactions. It 
is pointed out in [49J that one needs a sizable penguin weak annihilation amplitude in order to 
account for the data of B — > KqTt in the theoretical approach such as QCD factorization. 

Taking the cue from the constructive interference in the Cabibbo-allowed decay D + — > i^Q°7r + 
and noting that the phase of 02/01 in B — > Dtt decays lies in the first quadrant, one may be 
tempted to claim that the B~ — > Dq°7t + decay should have a rate smaller than the B° —¥ Dq~it + 
decay owing to a destructive interference in the former. This conjecture is supported in both the 



CLF model and heavy quark symmetry 78f] . The experimental observation that the production of 



broad D** states in charged B decays is more than a factor of five larger than that produced in 



neutral B decays (see Tables V and VI in [78J) is thus astonishing. This enigma as pointed out 



several years ago in 78] still remains unresolved. 



B(D+ -> if 1 (1270)vr+) is predicted to be 3.8 x 1CT 3 , 1.52% and 3.21% in [23J for 6 Kl = 33°, 45° and 57°, 
respectively. However, as noticed in Sec. III.B, the mixing angle 0K t has to be negative in the convention 
of QCD sum rules or the ISGW model. Therefore, the mixing angle chosen by [23[ does not have a correct 
sign. Also the relative sign between T and C amplitudes given there is erroneous. 
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VIII. CONCLUSIONS 



In this paper, we have studied the charmed meson decays into final states containing one 
pseudoscalar meson and one even-parity meson, the data of which are inferred from detailed Dalitz 
analyses of three-body decays and the finite width effects. The non-perturbative flavor diagram 
approach and the factorization calculation are undertaken to analyze these decay processes. In 
the diagrammatic framework, we have extracted the sizes and relative strong phases of various 
flavor diagrams in a least model-dependent way, based on current experimental measurements of 
decay rates. In the factorization approach, we have neglected weak annihilation diagrams {E and 
A) while making our predictions. Besides, we use the form factors evaluated in the covariant 
light-front model to compute the D decay branching fractions. 

In the D — > SP decays with S being in the non-strange nonet, our fits present several robust 
features against the 2- and 4-quark pictures. First, as in the PP and VP decays, the weak 
annihilation diagrams are non-negligible. However, it is a theoretically challenging puzzle that the 
A amplitude here is the largest one. Secondly, the relative strong phases 8a and 5e are preferably 
around 30° and 160°, respectively. Finally, the data do not prefer either the 2-quark scheme or the 
4-quark scheme. This conclusion can be best seen from a comparison between the D°' + — > fair and 
<77r decays, and also revealed in the Xmin vames i n Table IVIII1 We note that due to a paucity of 
measured modes, the diagrammatic approach is ineffective for the D — > KqP decays or the D — > AP 
decays. Theoretical calculations based on factorization have been applied to both D — > SP and 
AP transitions. The factorization hypothesis seems to work pretty well for those Cabibbo-favored 
modes that involve only the color-allowed tree (T) and color-suppressed tree (C) amplitudes. The 
factorization calculations for the other decay modes are typically about one order of magnitude 
smaller than experiment, conceivably due to the negligence of weak annihilation contributions. 

Contrary to the Cabibbo-allowed D + — > PP, VP decays where T and C amplitudes interfere 
destructively, the color-allowed and color-suppressed tree amplitudes in the Cabibbo-allowed decays 
D + — > K*° tt + , K^(1400)tt + contribute constructively. This explains why their branching fractions 
are large, of 0(2%). 

The D — > TP measurements poise the biggest problem for theory. Even though the magnitudes 
and phases of some amplitudes can be extracted from data, as given in Table IXIl quite opposite 
conclusions are reached when different sets of data [Fit (A) versus Fit (B)] are used. This could 
be caused by the D + — > K*2~k + and D + — > K^n^ decays, as explained in the text. Predicted 
branching fractions based on factorization are at least two orders of magnitude smaller than data, 
even when the decays are free of weak annihilation contributions. We cannot find possible sources 
of rate enhancement. 

We also examine the finite width effects for decays that are kinematically forbidden if the width 
of the even-parity meson is not taken into account. We find that the branching fraction of D + — > 
cnr + extracted from three-body decays is enhanced by a factor of 2, whereas B(D° — > j 2(1270) K°) 
is reduced by a factor of 4 by finite width effects. 

Our study shows that some of the above-mentioned puzzles call for measurements of yet ob- 
served decay modes as well as more precise determination in the decay rates. The other puzzles, 
particularly in the TP modes, demand better understanding of the underlying dynamics. 
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